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LECTURE 23. THE DIRAC DISTRIBUTION 



Impulse signals: Dirac's idea. Let a > be small. Let fa{t) be identically zero except for the 

/oo pa 
fa{i)dt = / fa{t)dt / 0. If the integral is not very small then fa{t) 
oo Jo 

must be quite large in the interval t € [0, a], and the function describes the "impulsive" behavior. 
An "impulsive function" means a signal which acts for a very short time but produces a large 
effect. The physical situation is exemplified by a lightning stroke on a transmission line or a 
hammer blow on a mechanical system. 



fait) = 




Figure 23.1. Graph of a typical unpulsive function 



In the early 1930's the Nobel Prize winning physicist P. A. M. Dirac developed a controver- 
sial method for dealing with impulsive functions. Let a 0+. The function fa{t)/{f fa{t)dt) 
approaches to, say, 6{t) which takes zero for t / and the integral of 5{t) over any interval con- 
taining is the unity. The function 5{t) is called "the Dirac delta function". 

To formulate the Dirac delta function, let a > and let 

1/a, te[0,a) 
0, elsewhere. 



fa{t) 



so that 
(23.1) 



fait) dt = 

* 6{t) as a 



fait) dt 

0+, then 



1 for all a > 0. 



Figure 23. 2. Graph of fa{t) 



6{t) = for t^O, 



6{t) dt = l. 



(This is often used as the definition of S{t) in elementary differential equations textbooks.) 
It is easy to show that 



nfait)] -- 

In this sense, S{t) describes the effect of 
(23.2) 



''-dt 



as a 



0. 



sa 



cm] = 1. 



Let us say that no ordinary function with the property (23.1) exists, and whatever else 6{t) may 
be, it is not a function of t\ However, says Dirac, one formally treats 6{t) as if it were a function 
and gets the right answer. 



1 



In the late 1940's, the French mathematician Laurent Schwartz* succeeded in placing the delta 
function on a firm mathematical foundation. He accomplished this by enlarging the class of all 
functions so as to include the delta function, called it the class of distributions. 

Here, we first explore the usefulness of d{t) in (23.1), and then we indicate its mathematical 
meaning. 

Examples. Let us consider the initial value problem 

y" + y = 6{t), y(0) = y'(0) = 0. 

Taking the transform, 

J^y = 2 ^1 1 ' y{t) = h{t) sin t. 
+ 1 

Here, we use (23.2). 

The solution y is continuous for alH G (—00,00) and it satisfies the differential equation every- 
where except t = 0. At i = 0, however, it satisfies neither the differential equation nor the initial 
condition. It is not even differentiable at t = 0. Indeed, y'{0+) = 1 and y'(0— ) = 0. The unit 
impulse signal 6{t) produces a jump of magnitude 1 in y'{t) at t = 0. 

Let us now consider 

y" + y = fM = .(0) = .'(0) = 0. 

10, elsewhere. 
We write fa{t) = -{h{t) — h{t — a)). Taking the transform then gives 



a 



+ 1 a 



and thus. 



Vaif) = -h{t){l — COS t) hit — a)(l — cos(t — a)) 

a a 

0, (-00,0] 
1 — cos t 

, tG(0,a) 

cos(t — a) — cos t 

, 0,00). 

a 

As a — > 0+, the second interval t G (0, a) vanishes, and the third interval t G [a, 00) tends to [0, 00) 
and the functional expression on this interval tends to sin t. In other words, as a ^ the solution 
ya{t) gives the same solution as was obtained with 5{t). We note that an uncritical use of 5{t) as 
in (23.1) and (23.2) gives the correct answer as obtained by a conventional passage to the limit. 
Moreover, the conventional method is much more difficult. It is where the usefulness of 5{t) lies. 

Next, let us consider 

y" + 2y' + 2y = 5{t), y(0) = y'(0) = 0. 

Taking the transform, 

= 2 I o 10 = 7 A2 I 1 ' = Kt)e-^ sin t. 

This example illustrates another feature of impulsive signals that it induces a lasting effect. 



'Laurent Schwartz wass awarded a Fields medal (a mathematical equivalent of the Nobel prize) in 1950 for his 
creation of the theory of distributions. 
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+ 2D + 2 



input 



Figure 23.3. The effect of impulsive signals. 



Theory of distributions. We conclude the lecture with a very brief description of the germ of 
Laurent Schwartz's brilliant idea. 

A function is characterized by giving its value at each t. A distribution 6{t) is characterized, 
not by its value at t, but by giving its value on a suitable class of functions (p, called test 
functions. Test functions are assumed to have derivatives of all orders and to vanish outside of a 
finite interval. 

Let us define 6 as 



(23.3) 



6{t)(P{t) dt = 0(0) 



for any test function (j). We note here that we cannot speak of the value of 5{t) at t. The only 

/OD 
(5(t)0(t) dt. The distribution 6{t) is never used alone, but only in com- 
-oo 

bination with functions. 

If d{t) were a function and if the integral in (23.3) were an ordinary integral, then by a change 
of variables. 



5{t - c)(j){t) dt 



6{t)(p{t + c) dt. 



By (23.3) the right side gives the value of 4>{t + c) at i = 0, that is, </<(c). Now we take it as the 
definition of the left side. Thus, 6{t — c) is defined as 

5{t - c)(p{t) dt = 0(c) 

J — oo 

for any test function 0. 

Similarly, if 6'{t) were a continuously differentiable function, then an integration by parts would 
yield 



J — ( 



6'{t)(l}{t) dt = - S{t)(P'{t) dt. 

J — oo 

By (23.3) the right side is — (/)'(0), and again we take it to define 6'{t) as 

6'{t)(l){t) dt = -0'(O) 

for any test function (p. With the same line of thought, 5^") {t — c) is defined as 

/oo 
(5(")(t-c)0(t)dt = (-l)>(")(c) 
-oo 

for any test function 0. 
For a < b, let us define 

/ 5{t-c)ct){t)dt= I {h{t-a)-h{t-b))5{t-c)(p{t)dt=Y^'^^' 

J a J —oo I 0) 



if a ^ c ^ 6 
otherwise. 



Lecture 23 



3 



18.034 Spring 2009 



With the choice (j){t) = e where s is a constant, we obtain 

e^'^^dit -c)dt= { 

[0, c<0 

when c ^ 0. When c < the result is zero. If c = then the above gives the formula C6{t) = 1, 
which agrees with (23.2). 

These are definitions, but the discussion shows that the definitions are consistent with the ordi- 
nary rules of calculus. This is the reason why 6{t) can be treated as a function of the real variable 
t even though it is not a function. 

Finally, assume that y{t) = for t < and y'{t) = 5{t) for t > 0. The Laplace transform suggests 
that sCy = 1, and in turn, y agrees with the Heaviside function h{t) except perhaps at i = 0. (At 
t = 0, neither the physics nor the mathematics of the problem is clear.) In this sense, 

h'{t) = 6{t). 



mt-c)]= / 

Jo 
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